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Abstract: — We present in this paper a compilation of the present status of experimen- 
tal delayed-neutron precursor data; i.e. /3-decay half-lives {T1/2) and neutron emission 
^ I probabilities (Pn) in the fission-product region (27 < Z < 57). These data are com- 

CN ■ pared to two model predictions of substantially different sophistication: (i) an update 

P5 , of the empirical Kratz-Herrmann formula (KHF), and (ii) a unified macroscopic- 

^ ^ I microscopic model within the quasi-particle random-phase approximation (QRPA). 

Both models are also used to calculate so far unknown T1/2 and values up to 
. Z = 63. A number of possible refinements in the microscopic calculations are sug- 

gested to further improve the nuclear-physics foundation of these data for reactor and 
astrophysical applications. 
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o 

(N ■ INTRODUCTION 

o ■ 

I Half-lives (T1/2) and delayed-neutron emission probabilities (Pn) are among the easiest 

measurable gross /3-decay properties of neutron-rich nuclei far from stability. They are not 
only of importance for reactor applications, but also in the context of studying nuclear- 
>^ [ structure features and astrophysical scenarios. Therefore, most of our recent experiments 

performed at international facilities such as CERN-ISOLDE, GANIL-LISE and GSI-FRS 
^ ■ were primarily motivated by our current work on r-process nucleosynthesis. However, it is 

^ , a pleasure for us to recognize that these data still today may be of interest for applications 

in reactor physics, a field which we practically left shortly after the "Specialists' Meeting 
on Delayed Neutrons" held at Birmingham in 1986. 



X 



I Our motivation to put together this new compilation of /3-decay half-lives and /3-delayed 

neutron-emission came from recent discussions with T.R. England and W.B. Wilson from 
LANL about our activities in compiling and steadily updating experimental delayed-neutron 
data as well as various theoretical model predictions (Pfeiffer et al, 2000). They pointed out 
to us, that their recent summation calculations of aggregate fission-product delayed-neutron 
production using basic nuclear data from the early 1990's (Brady, 1989; Brady and England, 
1989, Rudstam 1993) show, in general, that a greater fraction of delayed neutrons is emitted 
at earlier times following fission than measured. As a consequence, the reactor response to 
a reference reactivity change is enhanced compared to that calculated with pulse functions 
derived from measurements (Wilson and England, 2000). Therefore, the use of updated 
and values is expected to improve the physics foundation of the basic input data used 
and to increase the accuracy of aggregate results obtained in summation calculations. 

Since the tabulation of Brady (1989) and Rudstam (1993), about 40 new P^ values have 
been measured in the fission-product region (27 < Z < 57), a number of delayed-neutron 

^E-mail address: Bernd. PfeifferOuni-mainz.de 
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Figure 1: Chart illustrating the data available in the fission-product region. The new data 
evaluation represents a significant extension of measured Pn values. Some data in the old 
data set are not present in the new data set. 

branching ratios have also been determined with higher precision, and a similar number of 
ground-state and isomer decay half-lives of new delayed-neutron precursors have been ob- 
tained. These data are contained in our compilation (Table and are compared with two 
of our model predictions: (i) an update of the empirical Kratz-Herrmann formula (KHF) 
for /3-delayed neutron emission probabilities P^ and /3-decay half-lives T1/2 (Kratz and 
Herrmann, 1973; Pfeiffer, 2000), and an improved version of the macroscopic- microscopic 
QRPA model (Moller and Randrup, 1990) which can be used to calculate a large number of 
nuclear properties consistently (Moller et al., 1997). These two models, with quite different 
nuclear-structure basis, are also used to predict so far unknown T1/2 and P^ values in the 
fission-product region (see Table |lj) . 

EXPERIMENTAL DATA 

Most of the new /3-decay half-lives of the very neutron-rich delayed-neutron precursor 
isotopes included in Table || have been determined from growth-and-decay curves of neu- 
trons detected with standard neutron-longcounter set-ups. As an example, the presently 



B. Pfeiffer, K.-L. Kratz, P. Mdller, Delayed Neutron-Emission Probabilities 



3 



Table 1: Experimental /?-decay half-lives Ti^2 /3-delayed neutron-emission probabilities 
Pn compared to three calculations. 
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t 



This experimental data point was added after the manuscript was completed and is therefore not taken 
into account elsewhere in figures and tables. 
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0.985 


538 


5 


















550 


10 








1.65 


18 








280 


30 


216 


146 


146 


2.2 
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Pn 
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22.0 
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7.02 
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36 


1.092 
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69.824 


2282 


10 


920 
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12 


7.645 


15.749 


99.988 
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40 
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10.850 


99.728 
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20 


521 


351 


2347 


30 


9 


14.190 


42.231 


100.000 
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Exp. 
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64.646 


7 


123 


158 


158 


20 


10 
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0.000 


25 


1051 


(303 


G03 


0.12 


(J 


0.062 
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28.109 
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345 


345 
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96 


96 






9.276 


20.237 
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142 
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59.075 
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369 


306 
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359 
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0.006 
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0.049 
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Exp. 
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0.803 


0.361 


0.361 


9600 
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9440 


26147 


26147 
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4800 




4442 


11207 


11207 


0.000 


0.000 
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6821 


6821 


0.000 


0.000 
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1748 


3580 


3580 


0.000 


0.000 


0.000 






1226 


2527 


2527 


0.001 


0.000 


0.000 






764 


701 


701 


0.066 


0.020 


0.020 






570 


624 


624 


0.288 


0.469 


0.469 


10600 


1000 


9218 


40430 


40430 


0.000 


0.000 


0.000 






5219 


23562 


23562 


0.000 


0.000 


0.000 






2844 


12047 


12047 


0.001 


0.000 


0.000 






1794 


7521 


7521 


0.117 


0.144 


0.144 



Pn (%) 

Exp. KHF QRPA-1 QRPA-2 



used Mainz 47r neutron detector consists of 64 '^He proportional counters arranged in three 
concentric rings in a large, well-shielded paraffin matrix (Bohmer, 1998) with a total effi- 
ciency of about 45 %. The majority of the new P^ values were deduced from the ratios of 
simultaneously measured /3- and delayed-neutron activities. It was only in a few cases that 
7-spectroscopic data were used to determine the one or other decay property (e.g. indepen- 
dent Pn determinations for ^^Br, ^'^'^Rb and ^'^''Sn). Most of the new data were obtained 
at the on-line mass-separator facility ISOLDE at CERN (see, e.g. Fedoseyev et al, 1995; 
Kratz et al, 2000; Hannawald et al, 2000; Koster, 2000; Shergur et al, 2000). Data in 
the Fe-group region were obtained at the fragment separators LISE at GANIL (Dorfler et 
al., 1996; Sorhn et al, 2000) and FRS at GSI (Ameil et al, 1998; Bcrnas et al, 1998), and 
at the LISOL separator at Louvain-la-Neuvc (Franchoo et al, 1998; Wcissman et al, 1999; 
Mueller et al., 2000). Data in the refractory-element region were measured at the ion-guide 
separator IGISOL at Jyvaskyla (Mehren et al, 1996; Wang et al, 1999). Finally, some new 
data in the ^^^Sn region came from the OSIRIS mass-separator group at Studsvik (Korgul 
et al., 2000; Mach et al., 2000). 

In a number of cases, "old" Pn values from the 1970's deduced from measured delayed- 
neutron yields and (questionable) fission yields not yet containing the later well established 
odd-even effects, were - as far as possible - corrected, as was also done by Rudstam in 
his 1993 compilation (Rudstam, 1993). In those cases, where later publications explicitly 
stated that the new data supersede earlier ones, the latter were no longer taken into account. 
Multiple determinations of the same isotopes performed with the same method at the same 
facility by the same authors (e.g. for Rb and Cs precursors) were treated differently from 
the common practice to calculate weighted averages of experimental values, when a later 
measurement was more reliable than earlier ones. Finally, a number of "questionable" Pn 
values, in particular those where no modern mass model would predict the {Qp - S^) window 
for neutron emission to be positive (e.g. ^^^'^"^^Ba and ^^^La), are still cited in our Table, 



B. Pfeiffer, K.-L. Kratz, P. Mdller, Delayed Neutron-Emission Probabilities 



14 



but should in fact be neglected in any application, hence also in reactor calculations. 

MODELS 

Theoretically, both integral /3-dccay quantities, T1/2 and Pn, are interrelated via their 
usual definition in terms of the so-called /3-strength function (Sp{E)) (see, e.g. Duke et al. 
(1970)). 

l/ri/2= Sp{Ei)x f{Z,Qp-Ei)- (1) 

Ei>Q 

where is the maximum /3-decay energy (or the isobaric mass difference) and f{Z, Qp—Ei) 
the Fermi function. With this definition, niay yield information on the average (5- 

feeding of a nucleus. However, since the low-energy part of its excitation spectrum is 
strongly weighted by the energy factor of /?-decay, / ~ ((5/3 — Ei)^ , P\i2 is dominated by the 
lowest-energy resonances in Sp{Ei); i.e. by the (near-) ground-state allowed (Gamow- Teller, 
GT) or first-forbidden (ff ) transitions. 

The /3-delayed neutron emission probability (Pn) is schematically given by 

^ ^ j:f^S^iE,)fiZ,Q^-E,) 

YS' Sp{E,)f{Z,Qp-E,) 

thus defining Pn as the ratio of the integral /3-strength to states above the neutron separation 
energy S^^. As done in nearly all P^ calculations, in the above equation, the ratio of the par- 
tial widths for 1-wave neutron emission iJ'i^{Ey^)) and the total width (Ftot = r^i(-E'n) + 1^7) 
is set equal to 1; i.e. possible 7-decay from neutron-unbound levels is neglected. As we 
will discuss later, this simplification is justified in most but not all delayed-neutron decay 
(precursor - emitter - final nucleus) systems. In any case, again because of the {Qp — E)^ 
dependence of the Fermi function, the physical significance of the Pn quantity is limited, 
too. It mainly reflects the /3-feeding to the energy region just beyond ^n. Taken together, 
however, the two gross decay properties, Tx/2 and Pn, may well provide some first informa- 
tion about the nuclear structure determining /3-decay. Generally speaking, for a given Qp 
value a short half-life usually correlates with a small Pn value, and vice versa. This is actu- 
ally more that a rule of thumb since it can be used to check the consistency of experimental 
numbers. Sometimes even global plots of double-ratios of experimental to theoretical Pn 
to Tij2 relations are used to show systematic trends (see, e.g. Tachibana et al. (1998)). 
Concerning the identification of special nucleax-structure features only from T-^j^, and Pn, 
there are several impressive examples in literature. Among them are: (i) the development of 
single-particle (SP) structures and related ground-state shape changes in the 50 < < 60 
region of the Sr isotopes (Kratz, 1984), (ii) the at that time totally unexpected prediction 
of collectivity of neutron-magic (N=28) ''^S situated two proton-pairs below the doubly- 
magic ^^Ca (Sorlin et al, 1993), and (iii) the very recent interpretation of the surprising 
decay properties of ^^^'^32Q(i j^g^ ^bove A/" = 82 (Kratz et al, 2000; Hannawald et al, 2000). 

Today, in studies of nuclear-structure features, even of gross properties such as the T1/2 
and Pn values considered here, a substantial number of different theoretical approaches are 
used. The signiflcance and sophistication of these models and their relation to each other 
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should, however, be clear before they are applied. Therefore, in the following we assign the 
nuclear models used to calculate the above two decay properties to different groups: 

1. Models where the physical quantity of interest is given by an expression such as a 
polynomial or an algebraic expression. 

Normally, the parameters are determined by adjustments to experimental data and 
describe only a single nuclear property. No nuclear wave functions are obtained in 
these models. Examples of theories of this type are purely empirical approaches that 
assume a specific shape of ^^(E) (either constant or proportional to level density), 
such as the Kratz-Hermann formula (Kratz and Herrmann, 1973) or the statistical 
"gross theory" of (5-decay (Takahashi, 1972; Takahashi et al., 1973). These models 
can be considered to be analogous to the liquid-drop model of nuclear masses, and are 
— again — appropriate for dealing with average properties of /3-decay, however taking 
into account the Ikeda sum-rule to quantitatively define the total strength. In both 
types of approaches, model-inherently no insight into the underlying single-particle 
(SP) structure is possible. 

2. Models that use an effective nuclear interaction and usually solve the microscopic 
quantum-mechanical Schrodinger or Dirac equation. 

The approaches that actually solve the Schrodinger equation provide nuclear wave 
functions which allow a variety of nuclear properties (e.g. ground-state shapes, level 
energies, spins and parities, transition rates, T1/2, Pxn, etc.) to be modeled within 
a single framework. Most theories of this type that are currently used in large-scale 
calculations, such as e.g. the FRDM-I-QRPA model used here (Moller et al, 1997) or 
the ETFSI-|-cQRPA approach (Aboussir et al, 1995; Borzov et al., 1996), in principle 
fall into two subgroups, depending on the type of microscopic interaction used. An- 
other aspect of these models is, whether they are restricted to spherical shapes, or to 
even-even isotopes, or whether they can describe all nuclear shapes and all types of 
nuclei: 

(a) SP approaches that use a simple central potential with additional residual in- 
teractions. The Schrodinger equation is solved in a SP approximation and ad- 
ditional two-body interactions are treated in the BCS, Lipkin-Nogami, or RPA 
approximations, for example. To obtain the nuclear potential energy as a func- 
tion of shape, one combines the SP model with a macroscopic model, which then 
leads to the macroscopic-microscopic model. Within this approach, the nuclear 
ground-state energy is calculated as a sum of a microscopic correction obtained 
from the SP levels by use of the Strutinsky method and a macroscopic energy. 

(b) Hartree-Fock-type models, in which the postulated effective interaction is of a 
two-body type. If the microscopic Schrodinger equation is solved then the wave 
functions obtained are antisymmetrized Slater determinants. In such models, 
it is possible to obtain the nuclear ground-state energy as E =< *o|-ff|*o >) 
otherwise the HF have many similarities to those in category 2a but have fewer 
parameters. 

In principle, models in group 2b are expected to be more accurate, because the wave func- 
tions and effective interactions can in principle be more realistic. However, two problems 
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still remain today: what effective interaction is sufficiently realistic to yield more accurate 
results, and what are the optimized parameter values for such a two-body interaction? 

Some models in category 2 have been overparameterized, which means that their micro- 
scopic origins have been lost and the results are just paramerizations of the experimental 
data. Examples of such models are the calculations of Hirsch et al. (1992, 1996) where 
the strength of the residual GT interaction has been fitted for each element (Z-number) in 
order to obtain optimum reproduction of known and Pn values in each isotopic chain. 

To conclude this section, let us emphasize that there is no "correct" model in nuclear 
physics. Any modeling of nuclear-structure properties involves approximations of the true 
forces and equations with the goal to obtain a formulation that can be solved in practice, 
but that "retains the essential features" of the true system under study, so that one can 
still learn something. What we mean by this, depends on the actual circumstances. It may 
well turn out that when proceeding from a simplistic, macroscopic approach to a more mi- 
croscopic model, the first overall result may be "worse" just in terms of agreement between 
calcujlated and measured data. However, the disagreements may now be understood more 
easily, and further nuclear-structure-based, realistic improvements will become possible. 

PREDICTION OF P,, AND T1/2 VALUES FROM KHF 

As outlined above, Kratz and Herrmann in 1972 (Kratz and Herrmann, 1973) applied the 
concept of the /3-strength function to the integral quantity of the delayed-neutron emission 
probability, and derived a simple phenomenological expression for Pn values, later commonly 
referred to as the " Kratz- Hermann Formula" 

Pn ^ a[{Qp - S^)/{Qp - C)]' [%] (3) 

where a and b are free parameters to be determined by a log-log fit, and C is the cut-off 
parameter (corresponding to the pairing-gap according to the even and odd character of 
the /3-decay daughter, i.e. the neutron-emitter nucleus). 

This KHF has been used in evaluations and in generation of data files (e.g. the ENDF/B 
versions) for nuclear applications up to present. The above free parameters a and b were 
from time to time redetermined (Mann et al., 1984; Mann, 1986; England et al, 1986) as 
more experimental data became available. These values are summarized in Table |2|. Using 
the present data set presented in Table |^, we now again obtain new a and b parameters 
from (i) a linear regression, and (ii) a weighted non- linear least-squares fit to about 110 
measured Pn values in the fission-product region. For the present fits, the mass excesses to 
calculate Q13 and Sn were taken from the compilation of Audi and Wapstra (1995), other- 
wise from the FRDM model predictions (Moller et al, 1995). The cut-off parameter C was 
calculated according to the expressions given by of Madland and Nix (1988). With the con- 
siderably larger database available today, apart from global fits of the whole 27 < Z < 57 
fission-product region, also separate fits of the light and heavy mass regions may for the 
first time be of some significance. The corresponding fits to the experimental Pn values in 
the different mass regions are shown in Figs. 1-3, and the resulting values of the quantities 
a and b are given in Table |3|. It is quite evident from both the Figures and the Tables, that 
the new fit parameters differ significantly from the earlier ones; however, no clear trend 
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Figure 2: Fits to the Kratz-Herrmann- Formula in the region of "hght" fission products. 
The measured Pn values (dots) are displayed as functions of the reduced energy window 
for delayed neutron emission. The dotted line is derived from a linear regression, whereas 
the full line is obtained by a weighted non-linear least-squares procedure. For the fit 
parameters, see Table ^. 

with the increasing number of experimental data over the years is visible. With respect to 
the present fits, one can state that - within the given uncertainties - parameter a does not 
change very much, neither as a function of mass region, nor between the linear regression 
and the non-linear least-squares fit. However, for the slope-parameter b there is a difference. 
Here, the least-squares fit consistently results in a somewhat steeper slope (by about one 
unit) than does the linear regression. 

Based on the new non- linear least-squares fit parameters, the KHF was used to predict 
so far unknown Pn values between 27 Co and 63 Eu in the relevant mass ranges for each iso- 
topic chain. These theoretical values are listed in Table |l[ 

In analogy with the Pn values, the /3-decay half- lives T1/2 are to be regarded as "gross" 
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Table 2: Parameters from fits to the Kratz-Herrmann-Formula from literature. The two 
sets from Kratz and Herrmann (1973) derive from different atomic mass evaluations. 



Reference 


Parameters 




a [%] b 


Kratz and Herrmann (1973) 


25. 2.1 ±0.2 


Kratz and Herrmann (1973) 


51. 3.6 ±0.3 


Mann (1984) 


123.4 4.34 


Mann (1986) 


54.0 ±31/-20 3.44 ±0.51 


England (1986) 


44.08 4.119 



properties. Therefore, one can assume that the statistical concepts underlying the Kratz- 
Herrmann- formula for P„ values can be applied for the description of Ti/2- 




10"^ 10° 

Figure 3: Fits to the Kratz-Herrmann-Formula in the region of "heavy" fission products. 
For an explanation of symbols, see Fig. ^. 
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Table 3: Parameters from fits to tlie Kratz-Herrmann-Formula in different mass regions. 
Tlie sequence corresponds to Figs. 1 to 3. 



Region 


Lin. regression Least-squares fit 




a [%] b r^ 


a [%] b red. 


29 < Z < 43 


88.23 4.11 0.81 


105.76 5.51 80.97 
±37.67 ±0.61 


47 < Z < 57 


84.35 3.89 0.86 


123.09 4.68 57.49 
±41.17 ±0.38 


29 < Z < 57 


85.16 3.99 0.83 


80.58 4.72 78.23 
±20.72 ±0.34 



Tlie half-lives are inversely proportional to the Fermi-function f{Z, E), which, in first order, 
is proportional to the fifth power of the reaction Q^-value: 

T,/2^1/f{Z,E)^Qf (4) 




(Qp-Sn)/(Qp-C) 

Figure 4: Fits to the Kratz-Herrmann-Formula for all fission products. For an explanation 
of symbols, see Fig. ^. 
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Table 4: Parameters from fits to T1/2 of neutron-rich nuclides. 



lin. regression least-squares fit 



a 


ms] b r^ 


a 


ms] 


b 


red. 


2.74E06 4.5 0.72 


7.07E05 


4.0 


1.1E04 






±5.33E05 


±0.4 





Tlierefore, in a log-log plot of versus one expects the data points to be scattered 
around a line with a slope of about -(1/5). 

Pfeiffer et al. (2000) suggested to fit the T112 of neutron-rich nuclides according to the 
following expression: 

Ti/2 ^ax{Qp- C)-' (5) 

where the cut-off parameter C is calculated according to the fit of Madland and Nix (1988), 
and the parameters a and b are listed in Table |^. 

The gross theory has, basically, the same functional dependence on the Q^-value, but 
underestimates the /3-strength to low- lying states, which results in too long half- lives. We 
here compensate for this deficiency by treating the coefficient o as a free parameter to be 
determined by a fitting procedure. The values obtained are listed in Table |l|. 

PREDICTION OF T1/2 AND VALUES FROM FRDM-QRPA 

The formalism we use to calculate Gamow- Teller (GT) /5-strength functions is 
fairly lengthy, since it involves adding pairing and Gamow- Teller residual interac- 
tions to the folded- Yukawa single-particle Hamiltonian and solving the resulting 
Schrodinger equation in the quasi-particle random-phase approximation (QRPA). Be- 
cause this model has been completely described in two previous papers (Krumlinde 
et al., 1984; MoUer et al., 1990), we refer to those two publications for a full model 
specification and for a definition of notation used. We restrict the discussion here to 
an overview of features that are particularly relevant to the results discussed in this 
paper. 

It is well known that wave functions and transition matrix elements are more af- 
fected by small perturbations to the Hamiltonian than are the eigenvalues. When 
transition rates are calculated it is therefore necessary to add residual interactions 
to the folded- Yukawa single-particle Hamiltonian in addition to the pairing interac- 
tion that is included in the mass model. Fortunately, the residual interaction may 
be restricted to a term specific to the particular type of decay considered. To ob- 
tain reasonably accurate half-lives it is also very important to include ground-state 
deformations. Originally the QRPA formalism was developed for and applied only 
to spherical nuclei (Hamamoto, 1965; Halbleib et al., 1967). The extension to de- 
formed nuclei, which is necessary in global calculations of /?-decay properties, was 
first described in 1984 (Krumlinde et al., 1984). 

To treat Gamow- Teller j3 decay we therefore add the Gamow- Teller force 



(6) 
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to the folded- Yukawa single-particle Hamiltonian, after pairing has already been in- 
corporated, with the standard choice xgt = 23 McV/A (Hamamoto, 1965; Halbleib 
et ai, 1967; Krumlinde et ai, 1984; MoUer et ai, 1990). Here /3^±= EiO-itf are the 
Gamow- Teller /3='=-transition operators. 

The process of P decay occurs from an initial ground state or excited state in 
a mother nucleus to a final state in the daughter nucleus. For /3~ decay, the final 
configuration is a nucleus in some excited state or its ground state, an electron (with 
energy Ee), and an anti-neutrino (with energy E^,). The decay rate Wfi to one nuclear 
state / is 

y^fi = ^f-s\MM'fiZ,R,eo) (7) 
ri Zn'^ 

where R is the nuclear radius and eo = Eq/itlqC^, with itlq the electron mass. Moreover, 
|Myjp is the nuclear matrix element, which is also the /5-strength function. The 
dimensionless constant F is defined by 

r^^J^)' (8) 

where g is the Gamow- Teller coupling constant. The quantity /(Z, i?, eo) has been 
extensively discussed and tabulated elsewhere (Preston, 1962; Gove and Martin, 1971; 
deShaht and Feshbach, 1974). 

For the special case in which the two-neutron separation energy 6*211 in the daughter 
nucleus is greater than Q/j, the energy released in ground-state to ground-state (3 
decay, the probability for /3-delayed one-neutron emission, in percent, is given by 

Pin - 100 '""^"^^ (9) 
0<Ef<Qi3 

where Ef — — Eq is the excitation energy in the daughter nucleus and 5'in is the 
one-neutron separation energy in the daughter nucleus. We assume that decays to 
energies above 6*111 always lead to delayed neutron emission. 

To obtain the half-life with respect to /3 decay one sums up the decay rates Wfi 
to the individual nuclear states in the allowed energy window. The half-life is then 
related to the total decay rate by 

T, = ^^ (10) 

0<Ef<Qi3 

The above equation may be rewritten as 

J_ 27r^ln2 1 ^ B 

^ moc^ F2 ^ \Mfi\'f{Z,R,eo) E \MM'f{Z,R,eo) ^ ' 

0<Ef<Qf} Q<Ef<Qf} 



with 
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Figure 5: Calculated ^-strength function for ^^Rb in our standard model (Moller et al, 
1997). However, the deformation is not taken from the standard ground-state mass and 
deformation calculation (Moller et al, 1995). Instead the ground-state shape is assumed 
spherical, in accordance with experimental evidence. The figure shows the sensitivity of 
the calculated Pn value to small details of the model. Since there is no strength below the 
neutron separation energy, the calculated /3-delayed neutron-emission probability is 100%. 
However it is clear from the figure that just a small decrease in the energy of the large peak 
just above the neutron binding energy would drastically change the calculated value. 



For the value of B corresponding to Gamow-Teller decay we use 

B = 4131 s (13) 

The energy released in ground-state to ground-state electron decay is given in 
terms of the atomic mass excess M{Z, N) or the total binding energy Eund{Z, N) by 

Qf^- ^M{Z,N)-M{Z + 1,N-1) (14) 

The above formulas apply to the (3~ decays that are of interest here. The decay Q 
values and neutron separation energies S^^l are obtained from our FRDM mass model 
when experimental data arc unavailable (Moller et al, 1995). The matrix elements 
Mfi are obtained from our QRPA model. More details are provided elsewhere (Moller 
et ai, 1990). 

We present here two calculations, QRPA-1 and QRPA-2 of T1/2 and Pn- They 
are based on our standard QRPA model described above, but with the following 
enhancements: 
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Figure 6: This calculation corresponds to the QRPA-1 model specification. However, this 
nucleus is known to be spherical although a deformed shape was obtained in the ground- 
state mass-and-deformation calculation (Moller et al., 1995). Therefore, in our QRPA-2 
calculation in Fig. 0, this nucleus is treated as spherical in accordance with experiment. 

For QRPA-1: 

1. To calculate /3-decay Q- values and neutron separation energies Si,^ we 
use experimental ground-state masses where available, otherwise calculated 
masses (Moller et al, 1995). In our previous recent calculations we used the 
1989 mass evaluation (Audi 1989); here we use the 1995 mass evaluation 
(Audi et al, 1995). 

2. It is known that at higher excitation energies additional residual interac- 
tions result in a spreading of the transition strength. In our 1997 calculation 
each transition goes to a precise, well-specified energy in the daughter nu- 
cleus. This can result in very large changes in the calculated values 
for minute changes in, for example Si^, depending on whether an intense, 
sharp transition is located just below or just above the neutron separation 
energy (Moller et al, 1990). To remove this unphysical feature we intro- 
duce an empirical spreading width that sets in above 2 MeV. Specifically, 
each transition strength "spike" above 2 MeV is transformed to a Gaussian 
of width 

8-62 

= ^ (15) 

This choice is equal to the error in the mass model. Thus, it accounts 
approximately for the uncertainty in calculated neutron separation energies 
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Figure 7: This calculation corresponds to the QRPA-2 model specification. The calculation 
is identical to the calculation in Fig. ^ except that the ground-state shape here is spherical. 

and at the same time it roughly corresponds to the observed spreading of 
transition strengths in the energy range 2-10 MeV, which is the range of 
interest here. 

For QRPA-2: 

1. In this calculation we retain all of the features of the QRPA-1 calcula- 
tion and in addition account more accurately for the ground-state deforma- 
tions which affect the energy levels and wave-functions that are obtained in 
the single-particle model. The ground-state deformations calculated in the 
FRDM mass model (Moller et ai, 1992), generally agree with experimental 
observations, but in transition regions between spherical and deformed nu- 
clei discrepancies do occur. In the QRPA-2 calculation we therefore replace 
calculated deformations with spherical shape, when experimental data so 
indicate. This has been done for the following nuclei: 

67-78pg^ 67-79CO, 73-80^1^ 73-81^^, 78-842^, 79-S7q^^ 83-90Qg^ S^-gi^g^ 87-94gg^ 



87-96 



Br, 



92-98 



Kr, 



91-96 



Rb, 



96-97 



Sr, 



96-98 



Y, 



134-140 



Sb, 



136-141 



Te, 



137-142 



I, 



i4i-i43xe, and 



141-145 



Cs. 



To illustrate some typical features of /5-strength functions we present the strength 
function of ^^Rb calculated in three different ways in Figs. ^-|^. 

It is not our aim here to make a detailed analysis of each individual nucleus, 
but instead to present an overview of the model performance in a calculation of 
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Figure 8: Ratio of calculated to experimental /3-decay half-lives for nuclei in the fission- 
product region in three different models. 

a large number of /3-decay half-lives. In Figs. || and ^ we compare measured and 
calculated /5-decay half-lives and /9-delayed neutron emission probabilities for the 
nuclei considered here. To address the reliability in various regions of nuclei and versus 
distance from stability, we present the ratios T/3^caic/^/3,exp -Pn,caic/-Pn,exp versus the 
quantity Tg^exp • Because the relative error in the calculated half-lives is more sensitive 
to small shifts in the positions of the calculated single-particle levels for decays with 
small energy releases, where long half-lives are expected, one can anticipate that 
half-life calculations are more reliable far from stability than close to /3-stable nuclei. 

Before we make a quantitative analysis of the agreement between calculated and 
experimental half-lives we briefly discuss what conclusions can be drawn from a simple 
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Figure 9: Ratio of calculated to experimental /3-delayed neutron-emission probabilities for 
nuclei in the fission-product region in three different models. 

visual inspection of Figs. H and ||. As functions of Tg,exp one would expect the average 
error to increase as 7)3,exp increases. This is indeed the case. In addition one is left 
with the impression that the errors in our calculation are fairly large. However, this 
is partly a fallacy, since for small errors there are many more points than for large 
errors. This is not clearly seen in the figures, since for small errors many points are 
superimposed on one another. To obtain a more exact understanding of the error in 
the calculation we therefore perform a more detailed analysis. 

One often analyzes the error in a calculation by studying a root-mean-square 
deviation, which in this case would be 

C"rms = —^^{Ti3,exp — Tf^,calc) (16) 
^ i=l 
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Table 5: Analysis of the discrepancy between calculated and measured P -decay half-lives 
shown in Fig. M. 



Model 


n 


Mr, 






1 n 

a 

n 


/3,exp 

(s) 


KHF 


115 


-0.11 


0.77 


0.33 


2.15 


1 


QRPA-1 


115 


-0.02 


0.95 


0.50 


3.14 


1 


QRPA-2 


115 


0.13 


1.37 


0.61 


4.04 


1 


KHF 


187 


-0.40 


0.58 


0.41 


2.56 


all 


QRPA-1 


187 


-0.06 


0.87 


0.59 


3.88 


all 


QRPA-2 


187 


0.22 


1.67 


0.75 


5.75 


all 



However, such an error analysis is unsuitable here, for two reasons. First, the quan- 
tities studied vary by many orders of magnitude. Second, the calculated and mea- 
sured quantities may differ by orders of magnitude. We therefore study the quantity 
log(T^,caic/7>,exp), which is plotted in Fig. |, instead of (T^,exp - ^/j^cak)^- We present 
the formalism here for the half-life, but the formalism is also used to study the error 
of our calculated values. 

To facilitate the interpretation of the error plots we consider two hypothetical 
cases. As the first example, suppose that all the points were grouped on the hne 
T'/3,caic/T'/3,exp = 10. It is immediately clear that an error of this type could be entirely 
removed by introducing a renormalization factor, which is a common practice in the 
calculation of /?-decay half-lives. We shall see below that in our model the half-lives 
corresponding to our calculated strength functions have about zero average deviation 
from the calculated half-lives, so no renormalization factor is necessary. 

In another extreme, suppose half the points were located on the line T^,caic/^/3,exp = 

Table 6: Analysis of the discrepancy between calculated and measured /3-delayed neutron- 
emission probabilities P^ values shown in Fig. ^. 



Model 


n 


Mr, 


r'l 




< 


pmin 
n,exp 

(%) 


KHF 


86 


-0.31 


0.49 


0.36 


2.31 


1 


QRPA-1 


86 


-0.12 


0.76 


0.60 


4.02 


1 


QRPA-2 


86 


0.12 


1.34 


0.65 


4.51 


1 


KHF 


118 


-0.29 


0.51 


0.44 


2.76 


all 


QRPA-1 


118 


-0.18 


0.66 


0.62 


4.14 


all 


QRPA-2 


118 


0.11 


1.28 


0.75 


5.62 


all 
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10 and the other half on the hne T^,caic/^/3,exp = 0.1. In this case the average of 
log(T/3^caic/^/3,exp) would be zero. We are therefore led to the conclusion that there 
are two types of errors that are of interest to study, namely the average position of 
the points in Fig. ||, which is just the average of the quantity log(T/3_caic/^/3,exp), and 
the spread of the points around this average. To analyze the error along these ideas, 
we introduce the quantities 



r 




ri 


= logio('") 


Mr, 


1 " 
1=1 


ri 





0", 



10 



exp 



1 " 



^ i=i 



1/2 



(17) 



where is the average position of the points and is the spread around this av- 
erage. The spread can be expected to be related to uncertainties in the positions 
of the levels in the underlying single-particle model. The use of a logarithm in the 
definition of r\ implies that these two quantities correspond directly to distances as 
seen by the eye in Figs, ^j-^, in units where one order of magnitude is 1. After the 
error analysis has been carried out we want to discuss its result in terms like "on the 
average the calculated half-lives are 'a factor of two' too long." To be able to do this 
we must convert back from the logarithmic scale. Thus, we realize that the quantities 
Ml^ and a]^ are conversions back to "factor of" units of the quantities Mr, and 0"^, 
which are expressed in distance or logarithmic units. 



DISCUSSION AND SUMMARY 

In Tables ^ and ^ we show the results of an evaluation of the quantities in Eq. (P^) 
for Ti/2 and corresponding to [3 decay of the nuclei in table |l|. In the QRPA 
calculations the ratio between calculated and measured decay half-lives is close to 1.0. 
This shows, as pointed out earlier (Moller and Randrup, 1990) that no renormalization 
of the calculated strength is necessary. The mean deviation between calculated and 
experimental half-lives is a factor of 2-5 depending on model and half-life cutoff. Also 
the calculated Pn values agree on the average with the experimental data. Here the 
mean deviation between calculated and experimental data is a factor of 3-6, again 
depending on model and half-life cutoff. All half-life calculations agree better with 
data for shorter half-lives, cf. Fig. and Table |^. Therefore one can expect the 
models to perform better far from stability than what is indicated by the table. The 
/3-delayed neutron emission rates are also better calculated in the region of short half- 
lives and high Pn values, cf. Fig. | and Table |^. Again, this suggests calculated Pn 
values are more reliable far from stability than indicated by Table ^. 



B. Pfeiffer, K.-L. Kratz, P. Mdller, Delayed Neutron-Emission Probabilities 



29 



The KHF results appear more reliable than the QRPA results. This may seem 
surprising at first, because the KHF has minimal microscopic content compared to the 
QRPA. However, an advantage of the QRPA is that it provides so much more detail 
about /3-decay than does the KHF, namely the ft values of the individual decays, and 
the transition energies associated with those decays. A very detailed discussion of the 
possible sources of discrepancies between our QRPA results and experimental data is 
presented in Ref. (Moller and Randrup, 1990). One difficulty the calculations face is 
that the calculated half-lives depend on the energy of the transitions as (Q/3 — E)^. 
As an example we note that calculated half-lives for ^^Rb, for which — 9.28 MeV, 
change by a factor 1.5 for a change in transition energies by only 0.4 MeV. It is very 
difficult to reproduce transition energies to this accuracy in a global nuclear-structure 
model. 

For the QRPA-2 calculation we observe that the average of T^,caic/2]3,exp is consid- 
erably larger than 1, which corresponds to a correct average. One would have a priori 
assumed that this calculation would be in better agreement with experiment since 
we substitute calculated deformations for spherical deformations when so indicated 
by experimental data. However, since we do not include /3-strength due to forbidden 
transitions in our model, one would indeed expect that calculated half-lives be too 
low on the average. The non-spherical deformations that occur, contrary to experi- 
mental observations, in the QRPA-1 calculations in some sense simulate the missing 
low-lying forbidden /3-strength. However, a much more satisfying description would 
be to use correct ground-state deformations and develop some model to account for 
the strength related to forbidden transitions. 

The Pn values calculated in the QRPA-1 are on the average too low. At present 
we have no clear explanation for this result. An obvious correction to the model is 
to take competition with 7 emission into account, in particular for emission of /„ > 3 
neutrons. However, such a correction would further lower the ratio Tn,caic/^n,exp- 
One may speculate that an accounting for both this effect and forbidden transition 
strength in QRPA-2 would bring about satisfactory agreement. This possibility need 
to be investigated. 

We feel strongly that in a global, unified nuclear-structure model a single set of 
constants must be used over the entire chart of the nuclides, otherwise the basic 
foundation of the model is violated. However, for the purpose of generating the best 
possible data bases of half-lives and /5-delayed neutron-emission probabilities a com- 
plementary approach is reasonable. Just as we feel it is appropriate to use experimen- 
tal ground-state deformations, experimental single-particle levels, when known, could 
also be used as the starting point for the QRPA calculations. In practice the situation 
would be that in some regions, such as near the doubly magic ^^^Sn, many half-lives 
and Pn values would be unknown, but considerable information on single-particle level 
order and energies would be available. This experimental information could then be 
taken into account by locally adjusting the single-particle model proton and neutron 
spin-orbit strengths and the diffuseness of the single-particle well to obtain optimum 
agreement with the observed single-particle data such as the observed neutron single- 
particle sequence /7/2, P3/2, P1/2, and hQ/2 near ^^^Sn. The hope would be that the 
local agreement would be retained in some limited extrapolation away from the known 
region. Such a fairly limited extrapolation would be all that is required to reach the 
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isotopes in the fission-product region where experimental data are not yet available, 
cf. Fig. |I]. Limited studies along these lines have been undertaken by, for example, 
Hannawald et al. (2000). Other highly desirable enhancements to the calculations 
would be to include first-forbidden strength, perhaps first in a gross-theory approach 
and later from a new microscopic model. The cut-off parameter C in the KHF for- 
mula could be taken from the Lipkin-Nogami microscopic calculation instead of from 
the Madland-Nix macroscopic expression. The energy window {Qp — S^) could be 
reduced by 150 to max 500 keV to account for the angular-momentum barrier for 
emission of Z > 3 neutrons in for example ^^''I. 

In conclusion we note that we now have available about 40 new experimental T1/2 
and Pn values in the fission-product region. Data for additional nuclei in this region 
that are required as input in reactor criticality, astrophysical and other applications 
are provided from theoretical calculations. The substantial increase in available exper- 
imental data since the compilations by Brady (1989) and Rudstam (1993) is expected 
to have a significant impact on applied calculations. 
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